An Extension of Operational Calculus 

By JOHN R. CARSON 

THE Heaviside operational calculus postulates at the outset that 
the initial (boundary) conditions at reference time / = are 
those of equilibrium ; that is to say, the system is at rest when suddenly 
energized at time / = by a "unit" impressed force. By unit im- 
pressed force is to be understood a force which is zero before, unity 
after, time t = 0. 

In a paper published in Volume 7, 1929, of the Philosophical Maga- 
zine, Van der Pol briefly indicated the appropriate procedure for ex- 
tending the operational calculus to cover arbitrary initial conditions. 
The present paper is an exposition of this generalization for a system 
of a finite number of degrees of freedom, followed by an application to 
the differential equations of the transmission line. While stated in 
the language of electric circuit theory, it is to be understood that the 
processes are generally applicable to a wide variety of problems. 

We start with the canonical equations for a network of n degrees of 
freedom 

Zll-Jl + Z 12-^2 + ' • • + Zl nln = E\ 

(1) 

ZnJ\ +Z n Jl + • • • + Z n J n = E n 

where 



*-(^a+*+££*) 



(2) 



Now multiply the equations (1) by e~ pt throughout and integrate 
from to infinity. Also let J m and F m denote the Laplace transforms 
of I m and E m ; thus 

I m e~ vl dt, 

(3) 
F m = r E m e->» dt. 

Now let I m ° and Q m ° denote the initial values (at time / = 0) of 
I m and the charge Q m in the wth mesh; also let us replace z jk of (2) by 

z ih = pL ik 4- Rjk + i/pCik. (4) 
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We then have, replacing (1), the system of algebraic equations: 

Z\iJ\ + Z12/2 + • • • + «i»/« = F t + G\, 

(5) 

ZnlJl + z n2 J 2 + • • • + z nn J n = F n + G n . 

In these equations, G lt G 2 , • • •, G n denote the following summations: 

G x = Lull 9 - QiVCnP + ^12/2° - Qi°/C 12 p + • • • 

+ L ln I H ° - QnVClnp, 

(6) 

G n = L nl h - Qf/Cmp + • • • + L„„/„° - Q n °/C nn p. 

The right hand sides of equations (5) are thus known in terms of the 
impressed forces and the specified initial values of the currents and 
charges. They can therefore be solved in the usual manner for 
Ji, • ' ' f Jn- Thus 

_ F x + G, , F 2 + G 2 , ,_ F n + G n 

J m — 7 ' 7 I ' ' ' "1 y ' (i ) 

Having thus determined J\, •••,/„ as functions of p, I x , • • • , I n are 
determined as functions of time by the Laplace integral equation : 

MP) = f Im(t)erPt dt, p R > c, (8) 

Jo 

which completes the formal solution of the problem. Note that if 
Gi = G 2 = • • • = G n = 0, the solution reduces to the usual form. 

Equations (6) for G lf • • • , G„, may be written in a compact and 
elegant form as follows: Let 

T = hZHLiJih, 

u-IEZ-f-QO* (9) 

T is then the kinetic or magnetic energy stored in the network and U 
is the corresponding potential or electric energy. Then 

^-(SL-MSiL— ^ -•"■ (l0) 

The foregoing solution is compact, elegant and formally complete. 
In practical applications to networks of many degrees of freedom it 
may well present formidable difficulties in computation and interpre- 
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tation. This, however, is merely answerable to the complexity of the 
physical problem, and no simpler general solution can possibly exist. 

The foregoing method when applied to the differential equations of 
the transmission line, leads to the following differential equations 

(Lp + R)J = -±*+LP, 

» (11) 

(Cp + G)* = -±J+CV<>. 

Here / and $ are Laplace transforms of the current J and voltage V 
and 1°, V° are the initial values of I and V at reference time / = 0. 
/, 4>, 1°, V° are functions of x but of course independent of /. 
The formal solution of equations (11) is as follows: write 

Lp + R = Z(p) = Z, 

Cp + G = Y(p) = Y, (12) 



Also 



Then 



IZY = y, W? = K. 



LI°-^£-V° = F(x) - F. 
Y ox 



J = e-y* I A + ^ C*dyF(y)ey» J 

~ eyx { B + Tk r d y F ^ e ' yu } • 



(13) 



4> =-*!-/ + -£ F°. (14) 

7 dx Y 

A and B are constants of integration determined by the relations 
between / and * at the physical terminals of the line. 



